Strong field ionization to multiple electronic states in water 
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High harmonic spectra show that laser-induced strong field ionization of water has a significant 
contribution from an inner- valence orbital. Our experiment uses the ratio of H2O and D2O high 
harmonic yields to isolate the characteristic nuclear motion of the molecular ionic states. The nuclear 
motion initiated via ionization of the highest occupied molecular orbital (HOMO) is small and is 
expected to lead to similar harmonic yields for the two isotopes. In contrast, ionization of the second 
least bound orbital (HOMO-1) exhibits itself via a strong bending motion which creates a significant 
isotope effect. We elaborate on this interpretation by simulating strong field ionization and high 
harmonic generation from the water isotopes using the time-dependent Schrodinger equation. We 
expect that this isotope marking scheme for probing excited ionic states in strong field processes 
can be generalized to other molecules. 



The superposition of several ionic states resulting 
from strong field ionization of molecules has been 
observed via high harmonic generation (HHG). Tunnel 
ionization preferentially ionizes the HOMO. This leaves 
the molecule in the ground ionic state. Rotational 
or vibrational pre-excitation is typically used to force 
ionization from inner valence orbitals [TH3]. The 
population of ionic states is then probed by the 
recombination step of HHG 4 6 . In this contribution we 
show that strong field ionization of water can populate 
several ionic states without vibrational or rotational pre- 
excitation. This is done with an isotropic, unexcited 
ensemble by using isotope marked molecules. The 
excitation of lower lying valence states triggers an 
electronic wave packet with femtosecond to attosecond 
time scale dynamics 7 , which was demonstrated recently 
for the case of fine-split ionic states of Krypton [5] . 

We investigate ionization of the inner valence orbital 
by comparing the high harmonic yields from different 
isotopes of water. As first understood by Lein and 
observed by Baker et at, nuclear dynamics occurring 
between ionization and recombination are mapped onto 
the harmonic spectrum and can be isolated by using 
isotope marked molecules [TU]. Thus, the presence 
of multiple ionic states in HHG can be determined if the 
different ionic states launch significantly different nuclear 
motion. 

Figure [I] presents an intuitive picture of the theory 
elaborated in [9] . The ionization of the molecular orbital 
(either HOMO or HOMO-1 in Fig. Qa)) launches a 
vibrational wave packet on an ionic potential energy 
surface in Fig. 1(c) (red and green, respectively). The 
harmonic light is emitted as the returning electron 
recombines with the molecular ion. Photoionization and 
HHG recombination are inverse processes [TTj . The well 
known modulation of extreme ultraviolet photoionization 
cross sections by the Franck-Condon overlap (shown in 



Fig.Qb) for H 2 0) then requires the HHG recombination 
cross section to be modulated by the overlap of the ionic 
state nuclear wave packet and the neutral vibrational 
ground state. As the ionic state nuclear wave packet 
moves, it loses overlap with the ground state. The 
nuclear wave packet for the heavy isotope is slow, thus 
maintaining a stronger overlap than the lighter, faster 
isotope. This gives a larger harmonic yield for the heavier 
isotope. The timescale of the nuclear motion is mapped 
onto the energy of different harmonics [9: . Comparison of 
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FIG. 1: (a) Constant amplitude contours of the HOMO (l&i) 
and HOMO-l(3ai) orbitals. (b) Photoelectron spectrum of 
H2O reproduced from [12]. The transition to the ground 
ionic state ( 2 Bi) occurs near 12.6 eV and populates only 
a few vibrational lines. In contrast, the transition to the 
first excited ionic state ( 2 Ai), which occurs around 14.8 eV, 
populates many vibrational lines, (c) Potential energy curves 
of the neutral and first two ionic states. The minimum of 
the 2 Bi surface lies almost directly above that of the neutral, 
while the minimum of the 2 Ai surface is at 180°. 
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harmonic spectra from the different isotopes isolates the 
effect of nuclear motion, since their electronic structure 
is identical. 

In H2O, population of different ionic states can launch 
significantly different vibrational motion. Ionization 
of the HOMO, lb\, excites the molecule into the 
ionic ground state, X 2 Bi. The resulting potential 
energy surface [Fig. [ljc)] is very similar to that of the 
neutral as can be seen from the exponentially decaying 
Franck-Condon progression of the X 2 B 1 state in the 
photoelectron spectrum shown in Fig. [ijb). This is 
due to the non-bonding properties of the lbi orbital 
evident by its shape: the l&i is a single oxygen p orbital 
that sticks out of the molecular plane [Fig. [TJa)]. In 
turn, only extremely weak nuclear dynamics are induced 
[l"2l [T3] . In contrast, ionization of the HOMO-13ai, 
which sends the molecule into the A 2 Ai state, strongly 
excites the bending mode O |T3]. The equilibrium 
bond angle changes to 180° [FigJTJc)]. The 3ai has an 
intramolecular bonding character, see Fig. [TJa). The 
photoelectron spectrum illustrates the strong (weak) 
vibrational excitation upon ionization of the 3ai (l&i) 
orbital. The vertical binding energies of the 3ai and I61 
orbitals are 14.8 and 12.6 eV, respectively |14) . 

The very weak vibrational wave packet excited via 
population of the ionic ground state indicates that 
HHG from the I61 should not display a strong isotope 
effect. This is confirmed by a recent simulation, which 
considered ionization solely from the l&i and found 
a negligible difference between H2O and D2O high 
harmonics for vibrationally unexcited neutral molecules 

Several experiments have shown that inner valence 
orbitals can contribute significantly to high harmonic 
yields [TH51 ITS]. We expect the inner valence 3ai orbital 
to contribute to the high harmonics of water. The broad 
vibrational wave packet launched via ionization of the 
3ai orbital causes a rapid decrease in the overlap between 
the ionic and ground state nuclear wave functions. This 
results in a reduction of the HHG yield and significant 
difference between H 2 and D 2 harmonics. Water 
is therefore an ideal test molecule to determine inner 
valence contributions via nuclear motion. Harmonic 
generation in water has been demonstrated before in both 
liquid and gas phase samples [TTj [T5] . 

We generate harmonics with a commercial TkSapphire 
laser (pulse duration 30 fs, pulse energy 100 //J, central 
wavelength 800 nm, repetition rate 1kHz). A 40 cm focal 
length mirror focuses the beam to a full width at half 
maximum of about 60 /im. The harmonics between 20 
and 70 eV pass through an Al filter (thickness 100 nm) 
onto a toroidal flat-field grating, which disperses the 
radiation onto an extreme ultraviolet sensitive CCD 
camera. The gas cell is approximately 0.5 mm thick 
and has entrance/exit holes of approximate diameter 
0.2 mm. Water and heavy water reservoirs kept at 
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FIG. 2: (a) Heavy water to water harmonic yield ratios for 
several intensities. The error bars include the standard error 
of the mean of the individual datasets, but are dominated by 
a systematic effect characterized by taking two different sets 
of spectra at different times, (b) Simulated ratios. Markers 
represent simulations including both the l&i and 3ai orbitals 
for several intensities. Lines represent simulations that solely 
include the I61. See text for description of the theoretical 
error bars. 



room temperature are connected to the cell by a fine- 
dosing valve to keep the cell gas pressure of about 
3mbar (number density 10 cm -3 ). The laser intensity 
is changed by inserting pellicle beam-splitters into the 
beam. An iris placed just upstream of the grating rejects 
the long trajectory harmonics |19j . 

Figure [2]ja) shows the experimental ratios of heavy 
water to water spectra for several intensities spanning 
l-2xl0 14 W/cm 2 . Ratios calculated using odd harmonic 
peak height or area had negligible differences. The error 
bars are dominated by a systematic effect characterized 
by taking two different data sets of spectra at different 
times. The ratio increases monotonically with photon 
energy for all intensities. In the range 20-50 eV, the ratio 
spans values from 1 to 1.15. Also, changes in intensity 
did not measurably affect the ratio of any given harmonic, 
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effectively keeping the slope of the curve constant. 

These results show that the D2O/H2O harmonic ratio 
gets progressively larger for longer electron travel times. 
Thus strong field ionization launches significant and 
fast nuclear motion in water. As discussed above, the 
photoionization studies as well as the disagreement with 
the l&i-based simulation in [15] suggest that ionization 
of the 3ai is the dominant source of this motion. 

To model the experiment, we extend the theoretical 
approach used in |15j by taking into account multiple 
electronic states from different molecular orbitals. 
Furthermore, we solve the time-dependent Schrodinger 
equation under the single-active-electron approximation 
(SAE) with fixed nuclei separately for each orbital. Our 
approach has been used successfully for H 2 [20] as well 
as N2, O2, and CO2 [21] • The resulting spectrum is 
then multiplied by the nuclear autocorrelation function 
to include the nuclear motion. 

We now briefly outline this theoretical method, 
for details see references [30] [31]. Field- free Kohn- 
Sham orbitals are obtained using the LB94 exchange- 
correlation potential [22]. They are expressed in a 
multi-centered basis [23], which consists of a set of 
(typically atom-centered) local spheres (bases) defined 
as the product of spherical harmonics and a radial part 
expressed in _B-splines. Molecular symmetry (for water 
C21,) is accounted for by implicitly generating equivalent 
local basis sets at symmetry-equivalent positions. A 
large central sphere positioned at the charge center 
of the molecule defines an additional large B-spline 
basis. It overlaps with all atomic-centered spheres 
and serves for an improved description of the chemical 
bonding and the molecular electronic continuum. The 
latter is discretized with a density determined by the 
size of the central sphere. The radius of the central 
sphere is set to r^ ax = 121.5 a. u., the maximum 
number of angular momenta used is l m ax = 12 and the 
experimentally determined equilibrium geometry (Ron 
= 0.958i, <HOH = 104.5° [24) is applied. The time- 
dependent electronic wavefunction \^ a (t)} is represented 
as the sum of the Kohn-Sham orbitals multiplied by time- 
dependent coefficients. The index a refers to the orbital 
that is initially occupied. All Kohn-Sham orbitals are 
propagated in the velocity gauge during their interaction 
with a 10 cycle cos 2 laser pulse centered at 800 nm. The 
HHG spectrum for a single molecular orientation is then: 

/+00 
dte l ^ujd a (t)f{t)\ 2 (1) 
-00 

where d a {t) = ($ a (t)\p + A(t)\^ a (t)). The canonical 
momentum p and the vector potential A(t) are expressed 
in a.u. The damping function f(t) reduces d a (t) at the 
extreme times to minimize edge effects in the Fourier 
transform. 

Since the experimental sample is unaligned, we average 



simulated spectra from many different angles between 
the molecular axis and the laser polarization. This is 
performed for each orbital separately. Both coherent 
and incoherent averaging give essentially the same ratio 
D20/H 2 0. Averaging over only the three main axes of 
water (including the main ionization axes for the HOMO, 
HOMO-1 and HOMO-2) approximates the ratio obtained 
with about 80 orientations very well. Following Lein [9], 
the influence of nuclear dynamics is included only by the 
autocorrelation function C a (t) defined as 

c a (t) = (MoM(t)) (2) 

where \4> a ) denotes the vibrational ground state of 
H2O and \4>a) denotes the time-dependent nuclear wave 
packet launched on H2 + . Data for the autocorrelation 
functions is taken from Reutt et al |13j . Finally, we form 
the incoherent sum of the orientation-averaged HHG 
contributions S^o/i^o) ( w ) from the H OMO S lbl (u) 
and the HOMO-1 S 3ai {uj) (the HHG contribution of the 
HOMO-2 I&2 is found to be negligible) multiplied by the 
respective autocorrelation function 

s (d 2 o/h 2 o)^ )= \C^°l H *°\u)\ 2 S a {uj) (3) 

a— lbi ,3ai 

where we map the time of the autocorrelation function 
to energy by using equation (12) in 9J. This restricts 
the spectrum to energies below the cutoff. The ratio 
D2O/H2O is then obtained by S^°\uj) / S^°\lo) . 

Uncertainty in the theoretical ratios was determined 
by varying several parameters. First, three different 
damping functions were used. In each case the high 
harmonic spectra are additionally convolved with a 
Gaussian distribution function for four different variances 
chosen to maintain the basic peak structure of the 
spectra. Finally, the ratios using both peak heights and 
areas were computed for each case. 

Figure [2jb) shows simulations using the same 
intensities as the experiment. In agreement with |15j . 
the isotope effect calculated using only the lbi orbital 
(shown with lines) is very small in comparison to the 
measurement. The two orbital simulation (shown with 
markers) gives a much larger effect, confirming our 
qualitative estimate that the 3a± plays an important role. 
For example, at 45 eV, the ratio predicted by the two 
orbital simulation, approximately 1.10, is much closer to 
the measured value of 1.15 than that of our single orbital 
simulation, 1.04. 

The intensity dependence of the harmonic spectra is 
also telling. As stated above, the experimental slope 
had no discernible dependence on intensity. In contrast, 
increasing the intensity in the single orbital simulation 
decreases the isotope ratio for each harmonic, effectively 
reducing the slope of the curve. This is expected from 
the single orbital picture due to the one-to-one mapping 
between travel time and isotope ratio. In agreement with 
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FIG. 3: Ratio of the calculated total harmonic yields from the 
3ai and 16i orbitals of H2O as a function of intensity. Inset 
shows calculated D2O/H2O harmonic ratios of the 3ai and 
l&i orbitals separately for an intensity of 2 x W 14 W/cm 2 . 



the measurement, the two orbital simulation does not 
exhibit a discernible change in slope with intensity, even 
though the intensity range used is extremely large. 

Examination of the simulations suggests that the 
independence of the slope on intensity stems from an 
angular ionization saturation effect. Prior to the laser 
pulse, the angular probability distribution is isotropic. 
For a given molecule, the primary ionization direction 
of the lbi (perpendicular to the molecular plane) or the 
3ai (along the principle axis) are both equally likely to 
point along the polarization direction of the laser. As 
the intensity is increased, ionization causes the angular 
distribution to be depleted along the weakly bound I61 
ionization direction after several cycles of the pulse. The 
relative weight of HHG from the 3ai then increases with 
intensity (see Fig. [3]). The HHG spectrum of the 3<zi 
by itself provides an estimate of the upper bound for the 
isotope effect and is shown in the inset in Fig. [3] A 
maximum D2O/H2O ratio of approximately 1.4 results 
in this case. 

We have shown, using water, that strong field 
ionization populates an excited electronic state without 
vibrational/rotational pre-excitation. This is done with 
an isotropic ensemble by using isotope marked molecules. 
We modeled the multi-orbital system by employing a 
well tested SAE-TDSE simulation framework and show 
that both the lbi and 3ai orbitals contribute. The 
simulation predicts a magnitude of the isotope effect 
somewhat smaller than is measured but exhibits the same 
characteristic intensity dependence. The single orbital 
simulation predicts a much smaller effect and a different 
intensity dependence than is measured. The simulation 
suggests that the measured independence of the ratio on 
intensity is due to angular ionization saturation, which 



increases the relative contribution of the lower valence 
orbital at high intensities. 

These results provide a new tool for assessing 
contributions of different orbitals from strong field 
ionization. In addition, we show that ionization from 
inner valence orbitals can also excite nuclear dynamics. 
This extends the technique used in ^ [TU1 [25] to a 
larger range of molecules. Possible improvements of the 
theoretical model include, e.g., an incorporation of the 
relative phase between the two orbitals [2j and geometry- 
dependent ionization yields [26 . 
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